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ABSTRACT 


The small perturbation, two-dimensional transonic equation 
is manipulated with a separation-of-variables approach to 
obtain two ordinary, nonlinear, differential equations. 
Numerical integration of these differential equations results 
in new transonic boundary surfaces for planar external flows. 
A key ingredient in these solutions is the identification of 
dependence of two integrations constants, @ and £$, on the 
parameter (1-M,”). The anticipated behavior for both the Mach 
number and the pressure coefficient is used as a guide in the 
actual selection of the adjustable constants in the problem. 
The physical reality of our boundary surfaces is examined by 
displaying the boundary conditions they satisfy. The strictly 
sonic flow via O) has an analytic representation 
corresponding to a divergent surface which goes supersonic. 
This sonic solution is compared with an Euler-CFD approach 
confirming the validity of our results over the region where 
small perturbations apply. Solutions are also shown for 
M.=0.8,0.9,1.1, and 1.2 . These results are consistent with 
known behavior for both subsonic and supersonic external flow. 
Since the results of this work yield actual transonic 
contours, we can examine shockless surfaces for design 
applications. The possibility of starting with transonic 
surface is of interest to present day CFD approach. Finally an 
entire transonic upper surface is presented for M,=0.8, by 
patching a subsonic Mach number, which reaches a plateau at 
M=1.0, with a sonic flow. This patching requires the careful 
interpretation of a nondimensional reference length, called 


hoe woilch 1S a Lunction of M,. 
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a INTRODUCTION 

The governing nonlinear equations, together with realistic 
boundary conditions have always been the main difficulty in 
transonic aerodynamic analysis, [Ref. 1]. While the hodograph 
transformations approach linearizes these equations, the 
transformations are confined to a limited number of transonic 
flows due to complicated boundary conditions. The small 
perturbation approach simplifies the full potential 
Equations; furthermore, like other predictive methods, it 
succeeds in retaining the great sensitivity of transonic flows 
within small perturbations. This paper expands on the exact 
Selucion of the small perturbation, nonlinear, two 
dimensional, transonic equation, using the guidelines and 
tools given by Biblarz O., [Refs. 2 and 3]. Starting with a 
separation-of-variables approach, Chap. III reveals’ two 
nonlinear, ordinary, differential equations which lead into 
two exact solutions. Chaps. IV and V describe the numerical 
integration of these implicit solutions, which finally yield 
the boundary surfaces. Later, these boundary surfaces are 
found to satisfy the boundary conditions for a two dimensional 
surfaces. These findings are then compared with a finite 
difference solution of the Euler Equations for the sonic 
solution. The usefulness of the shockless, transonic, boundary 


surfaces appears once transformed to dimensional body 


surfaces; "patching" (1.e., translating), the dimensional, 
subsonic, body surface of C,=0 at the upstream inflow, with 
the dimensional, sonic body surface, enables us to form 


complete transonic upper surfaces of interest for design. 


Fre THE TRANSONIC EQUATION 


The small perturbation nonlinear, two dimensional 
Transonic Equation 
Z 
MM. +] 
G-m2)o, +o, = SeY*)) » ¢ (1) 


vy OE KX” KX 


is described by the velocity potentials 


u= ® v= @ (2) 


x M4 


Multiplying eqn. 1 by 1/U, enables us to transform the 


Transonic Equation to 
(1-M:) o,, + db, = Me (y+1)0,0,, (3) 


where the nondimensional velocity potentials are 


Oo, = b, - (4) 


me =a 
ua U,, 
The Meyer solution for the de-Laval nozzle defines a 


velocity potential 


z 2 4 
o(x,y) =b F~ + b(y+1) MA +b? (y+) ms (5) 


which satisfies the nondimensional form of equation (3) for 
M,=1.0 . The Meyer solution belongs to a general class of 


sonic flows as reported by Guderley and Yoshihara, [Ref. 4]. 


The transonic equation (1) may be modified to a more 
useful form if multiplied by M,?(y+1)/U.., instead of just 1/U- 
6, tO give 


(1-Mz) b,, + dy, = 0,0,,. (6) 


Here the modified velocity potentials are 


b, = Me (y+1) = d, = Me (y+1) (7) 
For a given problem, M, and y are constants, and eqn. 6 is 
really less cluttered. However, notice that M_7(y+1)/tUsmeee 
already given in the derivatives of @ in eqns. 7, and 
therefore appears in the expressions of C, and the boutdeme, 
surfaces aS will be shown later. This form of the modified 


transonic equation (6) will be examined in the next chapters. 


ely THE GOVERNING DIFFERENTIAL EQUATIONS 


An exact solution to the modified transonic equation 6 has 


been given by Biblarz, [Ref. 2], starting from 


b(x,yv) = oS(x,y) + (1-ME) x (8) 


where ¢°(x,y) refers to the sonic, small perturbation solution 
(M,j=1.0). In other words, "any flow which satisfies the sonic 
equation also satisfies the full transonic range with similar 
boundary conditions," [Ref. 2]. 

An exact solution may be obtained by using the separation- 


of-variables approach with the potential function ¢(x,y) of 


the form 
b(x,yv) = E(x)nly) + (1-MZ) x (9) 


Substituting the above @¢ function in the modified 
transonic equation 6 results in two ordinary, second order, 


nonlinear differential equations 








TE - 1 = 0 (10) 
= - An’ = 0 (11) 
y 


where ’} iS a separation constant. 


A simple solution to the first differential eqn. 10s 


obtained by multiplying both sides by dé/dx, 


Tela ele eee) ae _ (125 
(ze +e axe (*5) . 
OL 
Z/R) - Ae] 0 as 
Thus 3 
de _ | Bag: (14) 
ax 2 sage 
a s (15) 
34 52 + 
Z 
and 
dt 
Sa oe 
0 [- a (16) 
oo 


where @ and xX are integration constants. 
A similar procedure can be performed to solve the second 


differential eqn. 11 yielding 





in|) Go| an (4, 2) = Q (17) 
(<a dy *" 
Oms 
a fijdn\? _ Avs] . 18 
rabid 7 j — 


an = + 2h 4 + B (19) 


d 
an on 


(21) 
[Bees 


where $8 and y, are integration constants. 


Fay eel) 


Eqns. 16 and 21 are exact solutions to the differential 
eqns. 10 and 11. In fact, boundary conditions are needed to 
evaluate the integration constants a, 6, X, Yop aS well as A. 
An important step in this procedure is the recognition of a@ 


eam as functions of (1-M,”), [Ref. 2], i.e 


ae 


a = +C, (1-M2)° (22) 


B = +C, (1-MZ) (23) 


The upper positive sign and the lower negative sign in the 
above equations will be shown to belong to M,21.0 and M,<1.0 
respectively, namely, 

a20 & Ps0 imo}a Mesa 0 

a<0 & B<0 imoha Meade 


See Appendix B. 


The constants, a, 6 and A can be shown to be related by the 


expression «B ; 
Sad ee (aL 24 
4. ae ) (24) 
Since C, ,C, and } are positive constants 
_ 2 
CC, = =A (25) 


The usefulness of these definitions will be outlineqcmean 


more detail once we use the exact solutions from egns. 


CRP 


16 and 
while verifying the correct boundary conditions for the 


PECSSUNe COcr Meter ler 


Upon inserting eqns. 22 and 23 in eqns. 16 and 21 we have 


{- 3 = ———— 
en 22 (26) 
ici ac Gelinas 
a ee kc +f gy 
2X 3 2 (27) 
we ag Cal Vi 
Factoring out C, and C, we obtain 
_i dé 
ct = 3 eee SESE ee 
ee a Se [- (28) 





3A 2 BARR 
ac, ° + (1-MS) 


Introducing new dependent variables 








=i 
HI 
= 
—— 
NO 
~ 
peers 
wir 


the above simply becomes 


a aa 
ae. ‘ors (3) 2 if ae 
po eve 
1 ek 
_ 7 @ Alea S ay 
yor sh (S)* | 


Finally, if we define new independent variables 


(29) 


(30) 


(31) 


(32) 


(33) 


(34) 


(35) 


eqns. 32 and 33 then follow as 


dé 
xX = ee eee 
(36) 
Js E24 (1-m2)° 


an 


OO ———— 37) 
UD aes 5 


This set will be numerically integrated to obtain 
transonic surfaces with Cheir corresponding C, and Mach 1umbes 


profiles. 
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rV.. NUMERICAL INTEGRATION 


We can rewrite eqns. 36 and 37 in amore general form and 


specify the two integral limits 


2k 
x= f —— (38) 
Oza (heey - 
p AW 
Pea’ ge (39) 
ae BAS eae) 
¥(1-mMS) 3 


As indicated before, the upper and lower sign within the 
integrand root and the integral limits is used for M,21.0 and 
M.<1.0 respectively. In addition, our flow of interest is only 
the first quadrant in xX and Y, which represents physical 
external flows. To show this clearly, for the subsonic case 


MW<1.0, eqns. 38 and 39 reduce to the form 


xX = ee O>E>-(1-M.) (40) 


en, oe 


fj i 
0 i oi ee 2 (1-mMe) 3 (41) 
Ww? - (1-M2) 
(1-2) 3 


aids 


and for M,21.0 we have 


3 
0 


an, 


flz—(1-Mo) > (43) 


However, we will frequently retain the more general "+" 
notation in our discussions further on. 

Before we head to numerical solutions, we should indicate 
that sonic flow (M,=1.0) can be given in explicit form through 
eqns. 40,41,42,43, l1.e., the resulting equations for sonic 


flow are 





4 
AZ 
X = (44) 
{ 





.- <7 aw (45) 


3 
x (46) 


Ps Dae 
= (2 47 
ji (+ | (47) 


IZ 


Having defined £.4,X and Y in eqns. 30,31 and egns. 34,35, 


we can express the sonic solution rewriting eqns. 46 and 47 as 


1 32 
g = a) (48) 
22 
Ce a ie on 49 
(y-¥e)?( 25) yese 


Providing the sonic perturbation potential from eqn. 9 
° = En (50) 
and inserting eqns. 48 and 49 we have, [Ref. 2] 


me aes (51) 
2 ya 


Ss 


Next, eqns. 40 to 43 are numerically integrated and 
plotted in Figs. beans X29 fOr Sevarlous. Mes, (Esq; 
eee, 2. 0,1.1,1.2). Plotting 7 vs. Y for Y>0 in Fig. 2 1s an 
appropriate choice for investigating physical external flows 
as stated earlier. 

Our numerical solution has been based upon the Newton- 
Cotes method of order 4, [Ref. 5], which may be applied to an 


even number of subintervals, each described as, 


[eco dx = (22) 7 f(a) +32 £( 22) +12 £(252) +32 rte f(b) | 


(52) 


(See appendix A). 


a3 


Numerical solutions which express ¢ and n, aS has been 
noted, can be combined with egns. 30,31 and eqns. 34,35 and 
substituted in eqn. 9 to give a complete solution for the 
perturbation potential in terms of the positive constants A or 


(CG), welll ec ama) a 
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Vv. BOUNDARY SURFACES 


For an inviscid flow, the condition to be applied at the 
surface of a solid boundary is that the direction of the flow 
velocity vector is tangent to the solid surface, [Ref. 6]. In 
terms of perturbation velocities this boundary condition 


becomes 


(53) 


ie 
dx 


| V 
surface Us 


The modified velocity perturbation potential as pointed out 


earlier is 


ae: V 
>,, = M, (ytd) (54) 
thus 
fe we an by (55) 
oe M2 (y +1) 
We can show by substitution that 
(22) ee oe 
ax surface Mz (y +1) 


Now, deriving the perturbed potential function in eqn. 9 


with respect to y gives 


— ay 
Un =e a (57) 


gas. 


Eqn. 19 £Lor y=0 becomes 


on. | Bhp: (58) 
a i 


This relation can be rearranged using 4 and 6 in egns. 31 and 


23 to yield 


role 


ES 


~ 3 _ ue (59) 
ay i ie de) 


Upon uSing this equation and the relation of eqn. 30 to the 


derived velocity potential in egn. 57 results in 
Ws 
ob, = -§ 34) 2C,° ff? + (1-Mé) (60) 


The last equation can be rewritten with the benefit of the 


related constants C,, C, and A in eqn. 25, hence 
o, = -=2/ A? + (1-Me) (61) 


Once we have found the form of $,, Ehe expression forum. 


surface boundary condition is given as 


IL 


= of 
3 


ole 


) eee 


However, an important result to observe is that £ and 7 in the 
R.H.S. of this equation are given implicitly as a function of 
X and Y respectively in our earlier numerical integration 


solution, while the L.H.S. is expressed by x and y. To get 


is 


Sempatibility, we then transform the L.H.S. of eqn. 62 using 


th hain rule 
e chain ru dy _ a¥ dX/ax 


dx dx dyY/dy oo) 





Deriving eqn. 34 with respect to x and eqn. 35 with respect to 


y gives ; , 
oo : ce, (} (64) 
al i) 
Y 
therefore 


ho Se 2 a = 3 Me 
ay (+) PTT Ef 2 + (1-MS) (67) 
Oug 
dy B V2 i) 
SS oa (68) 
fi? + (1-MS) Meeaay, = 18) 


We are now in a position to compute the exact boundary 
surfaces out of eqn. 68. Starting with the L.H.S. we then 


define 


ey 


Having defined eqn. 69 we can plot E(Y¥) vs. Y in Fig. 3een 
various M,, (e.g., 0.8,0.9,1.0,1.1,1.2), Keeping in mindm@ems 
"+" notation for the supersonic (upper sign) and subsonic 
(lower Sign). 


The asymptotic value of E(Y) occurs at Y, where 


al 


A > #(1-M2) 3 7) 


In other words, the asymptotic value can be expressed as 
Yo=Y¥)(M.), or as function of (1-M,’), which has been 


numerically found to be 
a 


Y, = 5 
6 


(71) 
Ja -Me | 

Yo vs. M, is sketched in Fig. 4 to demonstrate an interesting 

symmetry with respect to M,=1.0 . Here Y=o but note that for 

practical reasons we can adopt finite values of Y, (perhaps as 

low as 7 to 20) to represent conditions sufficiently close to 

M,j=1.0 

Defining the integral of the R.H.S. of eqn. 68 as 


26 


2 
Ge = -{2| = (72) 
3/ M2 Pes 7") aes 


and redoing numerical integration allow us to plot K(X) vs. X 
for y=1.4 in Fig. 5. It is important to realize that Kies 


for M,<1.0 and K(X)<0O for M,21.0 , before solving eqn. 68. 


TEs: 


Therefore, integrating both sides of eqn. 68 yields for K(X) >0 


Yo 
[2H a¥ = K(x) (73) 
¥ 

and for K(X) <0 
N¢ 
fem d¥ = K(X) (74) 
Yo 


Once we specify X for K(X) we can immediately solve the L.H.S. 


integrals of eqns. 73 and 74 for Y, expressing 


(Y) a 1G) (75) 


surface 


namely, the required boundary surface for any chosen M,. 
We can rewrite eqns. 73 and 74 in the more complete form 


mote o<l.0 as 


4 xX 
{ dy i 1 fe ax (76) 
y ~ 2 S M: ( sf 

ame ror Mi=1.0 


Y - ; 6 

— = rT iH dx (72) 
vo yf? + (1-Mz) 31 Mo (y¥+1) 9 

Eqns. 76 and 77 are equivalent to eqns. 73 and 74 which 

finally yield the boundary surfaces. These are shown 

nondimensionaly in Figs. 6 and 7 for M,=0.8,0.9,0.97 and in 


ree 8 for Mo=1.0,1.1,1.2 


£9 


The sonic surface is most easily verified by direct 
analytical expansion of eqn. 77, using the sonic relations 
Stated in eqns. 46 and 47. 


Rie ae mei fi a oe 
{(3) oF -(4)atn [G42 - 


Yo 
Solving the above integrals on both sides gives 


yvi-y4i--_1 _/[2)\’x: (79) 
: (y+1)\3 


Recall that Y)y > » for the sonic case. However, choosing for 
example Y)=20 gives an asymptotic value of 7y=0.01, (eqn. 47), 
and E(¥,)=1000, (eqn. 69), for our explicit sonic solvtaam@ 
(See also Figs. 2 and 3). Therefore, dividing by a finite Yj 
enables us to obtain the nondimensional sonic boundary surface 
as a function of y 


| a (2)*(X) (80) 
(yl) \ Svea 


bs 


This will be shown to be a divergent surface which increases 


the Mach number to Supersonic conditions. 


20 


Vi. PRESSURE COEFFICIENT 


For a flow with small perturbations the pressure 


Se-erficient is given by, [Ref. 8] 


2 Z 
: $B (aa) 2 (#) 
U U Oe 


[2] +2) 


(81) 








A first order approximation for the linearized pressure 
coefficient in a two dimensional planar flow is 


Boe 
eae (82) 


Recall the modified velocity potential 


, = Me (y+2) — (83) 
Thus 
0, 
— ee See (84) 
U, Mo (yt) 
and from eqn. 82 we see that 
20. 
es “S555 SS (85) 
MeAy ts) 


Now, deriving the potential function in eqn. 9 with 


respect to xX we have 


>, = a + (1-MS) (86) 


Zell 


Rewriting eqn. 14 using the relation for a in eqn. 22 as 


5 
d 3A 
a = 5 es ri C, (1-M2)° (87) 


By taking advantage of £ as sitated ain "ec cniemee iemeeial factominre, 


out C,, eqn. 87 becomes 


. =o3 fa (-M)" (88) 


Substituting eqn. 88 and the relation for 7 in eqn. 32 Game 


velocity potential ¢, in eqn. 86 yields 


oe SE Cae) (89) 





Once again, using the benefit of the related constants C,, C, 


and A in eqn. 25, eqn. 89 becomes 
3 
o, = Hoy 2? + (1-M2)* + (1-m2) (90) 
therefore eqn > FOr se ee ele mas 


=? es 2 2,2 2 
ey Se) lee) (91) 
a aoe | 


We must anticipate the above expression for C, to Satisit aaa. 
boundary conditions on external boundary surfaces as will be 


CxplOored srhucener sone 


Vals. THE LOCAL MACH NUMBER 


The velocity vector V is given by the undisturbed uniform 


velocity U, and the perturbation velocities for the 2D flow as 
V=i(U.+u) + jv (92) 


Thus, we may obtain the local Mach number M in terms of the 


above velocities and the local speed of sound "a". 


eee Wile ee Ue (93) 


ee a? 


Neglecting higher powers of the perturbation velocities gives 


2 2U 
Ne cen 
. | u (94) 


oo 


Me = 
ae 


In addition, the ratio of the local speed of sound to the 


undisturbed uniform speed of sound becomes 





a* Th 
ge. a! ae (95) 
ae ‘ie 
Or 
2 ie) ie 
a = 0 (96) 
az cy patl 


were T, designates the stagnation temperature. 
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96 can be replaced using the well known relation 








Eqn. 
, 1+ tem 
=. = (97) 
ae Te ye 
2 
97 to egies 


Substituting the local speed of Sound tremgcenu: 
results as 
a+ 38) 


Niel =e 
(98) 


N 








98 is reduced to the form 


Having stated Cyan veqa se ecu. 
M? . Pe eee) 

—— y=. =e (99) 
ik M? 1a mM: 


Rearranging eqn. 99 and factoring out M yields the final 
expression for the local Mach number 


Me ) 
eee (100) 


M* = 
1+J21 Mc 
2 Pp 


Note that M=M, when C,=0 as it should. 
combined with the corresponding 


This important relation, 
will be used when discussing the boundary 


eqn. 94) One. 


conditions in the next chapter. 
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Ver BOUNDARY CONDITIONS 


Rather than providing the boundary conditions for the 
governing differential eqns. 10 and 11 to evaluate the 
constants, we imposed a solution by defining @ and 8 in egns. 
22 and 23 respectively. That is to say, boundary conditions 
have been already implied in the definitions of a and $. This 
Procedure led us to reveal the boundary surfaces upon using 
Mie tangency condation and compute the pressure coefficient C, 
on the body surface, (eqn. 91). 

Thus, we arrive at the necessity to verify the validity of 
aie eexpression for C, for the boundary surface, and verify the 
following requirements: 

® At the upstream end (X/Y,=0), we must require that C,=0 
ensuring the unperturbed Mach number of the transonic 

inflow to equal M,. 

® For the subsonic inflow at the downstream end (X @ Y/Y)=1) 
we must require that the local Mach number saturates 
towards M=1.0 . This is because of the well known relation 
in gas dynamics for a convergent flow. 


We first restrict our attention to the first requirement, 


MomeM,=1.0 at X/Y,=0 
E = e° = O (101) 


i 
ij oe ij ae Se ie 3 (102) 
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Put 


Inserting the above terms of € and % in egns. 90 and 91 gives 
as 
o, = -(1-M2) 3? \f(1-m2)* + (1-2) = 0 (103) 


and 
Gi @ 3/ Va8 0 es eee (104) 


Therefore, from eqn. 100 


M @. (eX? Vase oY) See ee (105) 


oo 


Hence, eqns. 104 and 105 comply with the first requirement. 
The complete behavior of the pressure coefficient and the 
local Mach number are plotted in Figs. 9 and 10 for M,=1.0 and 
M.=1.1. In addition, it 1s important to note that the boundary 
surfaces ought to be truncated at M,s1.2, so as not to exceed 
the small perturbation assumptions. 


Now, for M,<1.0 at X/Y,=0 


Fak" = 9 (106) 
However, at Y/Y,=0 
7-7 =” (107) 
Thus, eqn. 91 results as 
Cr @( X/Y, =O, Yvan = 203) Sone (108) 


We immediately conclude that the above coordinate (0,0) can 
not serve as the upstream start for the subsonic inflow. 
Instead, we must localize a new upstream end coordinate along 


~- 


the boundary surface whose £° and % values satisfy eqn. 91 


ZS 


melee equality C=O - Since every coordinate on the boundary 
sureace (X/¥,,Y/Y¥,) is a function of g and 9, (eqn. 76), we can 
numerically scan the locus of all the coordinates, i.e., the 
boundary surface, for such values of ~* and nH that will 
Satisfy eqn. 91. For example, scanning the results for the 
subsonic boundary surface of M,=0.8 reveals; £=£*=-0.35763 and 
n= =3.01393 at the coordinate (X/Y,=0.31366,Y/Y,=0.40616). 
Mis=rting the above values of — and % in eqn. 91 yields 
90003 . Hence, by relocating a new upstream coordinate for 
the subsonic inflow we are then able to satisfy our first 
requirement. (This is shown later in Figs. 11 and 12). 

We now consider the second requirement for the subsonic 


outflow at the downstream end where 


E- &* = -(1- Me) (109) 
a 110 
Gx y=) = = (lems)? ( ) 


or equivalently dY/dx > 0 


Substituting these terms for ~ and % in eqns. 90 and 91 gives 


ob, = (1-Mé) (laa) 

and ; 
c= D2 (ace) vet 

M2 (y+1) 
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Thus, we can express the local Mach number uSing eqn. 100 


ARMS) 
Mz (y +1) 
1- vest) M2 IES), 


2 ~ M2 (y+1) 


Z 


uals || A 


Meme (113) 


Upon expanding eqn. 113 and factoring out M,’” we obtain 


Me (y-1) +2 

oS eee (114) 
M2 (y-1) +2 

or equivalently 
M=1 (115) 


Hence, we have verified the second requirement: The local 
Mach number saturates at the downstream end to the value M=1.0 
for any subsonic inflow M,<1.0 . The pressure coefficient and 
the local Mach number behavior along the boundary surface are 
plotted in Figs. 11 and 12 for M,=0.8 and M,=0.9 respectively, 
Showing clearly both the upstream and downstream boundary 
Gonawtr onse 

Knowing that boundary conditions require the gradient of 
the potential function ¢ to vanish far away from the body 
gives rise to an important queStion: How is the flow pattern 
along the boundary surface propagated into the flow field? We 
treat this question by considering a computationally efficient 
Flow solver developed by J. A. Ekaterinaris, [Ref. 7]. This 
CFD solver consists of the inviscid Euler equations and uses 


fully implicit two dimensional Crank-Nicholson method with a 


Za 


central difference ADI scheme, and fourth and second order 
numerical dissipation, [Ref. 9]. This solver has been modified 
here to solve the inviscid flow around our earlier boundary 
Siierace SOlutiom for the chosen sonic case M,=1.0 . By 
creating a grid mesh which consists of a boundary surface 
edge, we are able to map the flow pattern. Upon observing the 
constant Mach lines as graphed in Fig. 14, we deduce that the 
flow is shockless through the whole transonic range; a shock 
is formed further downstream at the local Mach number M=1.9 

These results are in complete accordance with the inviscid 
small perturbation transonic solution, (see Fig. 9). In 
addition, each vector on Fig. 13 describes the resultant flow 
direction at a specified grid point, (not the streamlines !), 
colored by Mach levels. Thereby, we may evidence the far-field 


boundary conditions 


= 10 i = (116) 


In other words, the flow pattern 1S propagated out from the 
boundary surface gradually towards satisfying the far-field 
boundary conditions far away from the body. 

It is of special interest to inspect the dimensional 
boundary surface for possible design applications; For 
example, the dimensional, subsonic, boundary surface of 
M.=0.8 is obtained by multiplying the nondimensional 
coordinates (X/Yoo.s),¥/ Yoo.) by Yoog=2.84 ce — the 


dimensional, downstream end coordinate results as (X,,2.84), 
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as plotted in Fig. 15a. Similarly, multiplyinggied- 
nondimensional, sonic coordinates (X/Yoq9,¥/Youo) PY Youn=14, 
yields the dimensional, sonic, boundary surface, where the 
dimensional, upstream, starting coordinate becomes (0,14). 
Now, translating the dimensional, sonic, boundary surface, 
Upale( LSbae by (AX, AY) =(X,, - (14-2.84)), to match the 
downstream end coordinate of the above dimensional, subsonic, 
boundary surface, "patches" both dimensional boundary 
surfaces, to form a 2D planar upper surface as shown in Fig. 
15. We may use this technique to "patch" any subsonic boundary 
surface in order to form our transonic upper surface of 


interest. 


310) 


38S CONCLUSIONS 


We have shown that an exact solution for the transonic 
equation satisfies the complete behavior of the pressure 
coefficient, and the local Mach along a 2D-surface for planar, 
external flow. It is of special interest to inspect the 
dimensional boundary surface for possible design applications; 
For example, the dimensional, subsonic, and sonic boundary 
surfaces are obtained by multiplying the nondimensional 
coordinates (X/Yo,Y/Yo), by the corresponding Y,. Translation 
of the dimensional, sonic, boundary surface, to match the 
downstream end coordinate of the above dimensional, subsonic, 
boundary surface, "patches" both dimensional boundary 
surfaces, to form a 2D planar upper surface as shown in Fig. 
eS This technique may serve as a powerful tool for the 
implementation of supercritical, shockless, planar surface 
designs, because of the saturation of % with Y (Fig. 2), we 
can also describe certain internal flows such as planar 
converging, diverging tunnel walls. We may also extend these 


results to axisymetric transonic flow. 


a 
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APPENDIX A 
COMPUTER PROGRAMS 


CR KKK KKK KKK KKK KKK IKK KKK KKK KEK KEK KKK KKK KKK KKK KKK KEKEKKKKKKEKEK 


* THIS FORTRAN PROGRAM USES NUMERICAL INTEGRATION BASED ON 
THE NEWTON-COTES FORMULA TO SOLVE NUMERICALLY EQNS. 40 AND 42. 


THE OUTPUT IS IN THE FORM OF & vs. X 
KKK KKK KKK KKK KKK KKK KKK KKK KKK KK KKK K KKK KKK KKK KK KKK 
program KSI 
PARAMETER (N=5) 
IMPLICIT REAL*8 (A-H,0O-2Z) 
IMPLICIT REAL*8 (M) 
IMPLICIT REAL*8 (K) 
INTEGER INDEX, I 
DIMENSION X(N) ,XDOT(N) ,SAVEX(N) ,SAVED(N) , KSI (N) 
DIMENSION MACH (N) ,X1(N) , X2 (N) 


MACH (1) =0.8D0 

MACH (2) =0.9D0 

MACH (3) =1.0D0 

MACH (4) =1.1D0 

MACH (5) =1.2D0 

KSI (1) =0.0D0 

KSI (2) =0.0D0 

KS1I(3)=1.0D-3 

KSI (4) =0.0D0 

KSI (5) =0.0D0 
CORIO II IO II TO TO I TOI IOI FOI IO II IO IO IO a 40% 
Cc MAIN PROGRAM * 


CR KEK KKK KKK KKK KKK KKK KKK KKK KK KKK KKK KKK KKKKKK KKK KK KK KKK KKKKK 


DO 100 I=1,5 
he tsbr.3) THEN 
H=-0.001D0 
Gory s 
END IF 
H=0'. O1D0 
5 me tor .3) THEN 
TP=- (1.0D0-MACH (TI) **2) 


GOTO 10 
END IF 
ie — ih. ODO 


mee LF(L.LT.3.AND.KSI(1) .GT.TF.OR. 

1.GE.3.AND.KSI(I).LE.TF) THEN 
30 do 50 index=1,5 

CALL DERIVative(KSI,N,1I,X1,X2,XDOT, MACH) 

CALL NEWTONCOTES (KSI, XDOT,X, SAVEX, SAVED,H,N, INDEX, I) 
50 end do 
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WRETE (10425 ) 0) Se) 


GOTO 10 
END IF 
100 END DO 
= FOE 
END 
CORK IO KO tO OR Ok Kk 
G SUBROUTINE DERIVative * 


CORK KKK KKK KKK KKK KKK KKK KKK KEKE KKK KKK KKK KKK KKKKKKKKKKKKKKEKKEKEKE 


SUBROUTINE DERIVative (KSI,N,1I,X1,X2,XDOT, MACH) 
IMPLICIT REAL*8 (A-H,0O-2Z) 
IMPLICIT REAL*8 (M) 
IMPLICIT REAL*8 (K) 
DIMENSION XDOT(N) ,X1(N) ,X2(N) ,KSI(N) ,MACH(N) 
GO-TO, 420; 20,107 16710) ae 
10 X1(1)=KSE( aq) **2+ (12 0D0-MACH 442) 2 
K2 (1) =XiG) ss (a Obey Ss 40be) 
XDOl (i jai Ob, x2 am) 
GO he wee 


20 M1 (1I)=- (KSI(1) **2-(170D0-MAGH 1) ~4*2) 4-2) 
K2 (1) =- Xie O07 seeb a) 
ADOT CE) = 120072 an) 


CMOS Slchameneey) 

END 
CORR RO RR RR RR RK RR a a a RK KR KR aR KR KR aR KR KR KR KR KR KR KR KR KK KK KK KKK 
G SUBROUTINE NEWTONCOTES * 


COR RR KR KR KKK KKK KKK KKK KKK KEK KKK KEKE KKK KEK KKK KKK KKK KEKEKKKKEKKEKKEKEKE 


SUBROUTINE NEWTONCOTES (KSI, XDOT, X, SAVEX, SAVED, H,N, INDEX, I) 
IMPLICIT REAL*8 (A-H,0-2Z) 
IMPLICIT REAL*8 (M) 
IMPLICIT REAL*8 (K) 
DIMENSION X(N) ,XDOT(N) , SAVEX(N) , SAVED (N) , KSI (N) 
GO TO MG 2 3 7 4 5) eed enc 
i SAVEX (I) =X (TI) 
SAVED (1) =XDOT (T) 
KSI (I) =KSI (I) +0.25D0*H 
RETURN 
Z SAVED (I) =7.0D0*SAVED (I) +32.0D0*XDOT (T) 
KSI) KSI (1) +e. 25p04n0 
RETURN 
3 SAVED (2) =SAVED( 1) +22 700+ %DOln 
KS Lt E)=KSi (2) -0h eho 
RETURN 
4 SAVED (I) =SAVED (I) +32 .0D0*XDOT (T) 
KSI {(@)=KSi Tt) 202 25D645 
RETURN 
5 X (I) =SAVEX (I) + (H/90.0D0) * (SAVED (I) +7.0D0*XDOT (TL) ) 
RETURN 
END 


CORK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KEKE KEKKKKKKKK KK KKKKKKEKK 
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KKK KKK KKK KKK KKK KKK KKK IKK KKK KKK KKK KKK KKK KKKKKKK KKK KKK KKK KKK KKK 


THIS FORTRAN PROGRAM USES NUMERICAL INTEGRATION BASED ON 


THE NEWTON-COTES FORMULA TO SOLVE NUMERICALLY EQNS. 41 AND 43. 


THESOUTE UT IS stn tae FORM OF 7 vs. ¥ 


KKKKIKKKK KKK KKK KKK KEKE KKK KKK KE KKK KEKKKKKKKKKKKK KKK KK KKK KKKRKKKK 


program ETTA 

PARAMETER (N=5) 

PARAMETER (K=2000) 

IMPLICIT REAL*8 (A-H,0-2Z) 

IMPLICIT REAL*8 (M) 

Ni EGER INDEX, 1,7 

DIMENSION Y(N) ,YDOT(N) ,SAVEY(N) ,SAVED(N), 
PEVMENSTON BGAN), BWALI(K) vil (kK) ,Y1(N), Y2 
MACH (1) =0.8D0 

MACH (2) =0.9D0 

MACH (3) =1.0D0 

MACH (4)=1.1D0 

MACH (5)=1-2D0 

H=0.01D0 


(N) , MACH (N) 


TF 
(N) 


i] 


CK KKK KK KKKK KKK KKKKKKKKKKK KKK KKKKKKKKKKKKKKKK KK KKK KKKKKKR KK KKK 


& 


INTEGRAL FINAL LIMITS us 


CR KKK KKK KKK KK KKK KKK KKK KKK KKK KKKEKKKKKKKKKKKKKKKKKKKKKKKKKKKK 


i) 


BO / [=1,5 
PE MAGI) HO] LObO) -ToEN 
TF (I) =0.0D0 
Coton 
END IF 
TE (MACH (1) .=GT.120D0) > 2HEN 
TF(2)=(-(1.0D0-MACH (1) **2) ) **(1.0D0/3 .0D0) 
GOO a7 
END IF 
ia) (ae OO MACH (1) ** 2 es { 1 O03 700) 
END DO 


CKRK KK KKK KKK KKK KKK Kh KKK KKK KK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KK KK 


Se 


MAIN PROGRAM ‘ 


CHK KK KK KKK KK KKK KKK KKK KKK Kh KK KKK KKK KKK KKK KK KKK KKK KKK KKK KK KKKKK 


10 
20 


5/19) 


mee OO T=1_ 5 
Jem 
on (bP) =14..0D0 
mit) =2.0D0/7 (BTA(1) ** (1.0D0/2 .0D0)) 
ETF (ETA(1I)-H.GT.TF(1I)) THEN 
do 50 index=1,5 
CALL DERIVative (ETA,N,1I,Y1,Y2,YDOT, MACH) 
CALL NEWTONCOTES (ETA, YDOT, Y, SAVEY, SAVED,H,N, INDEX, 1) 
end do 
Bee Ggi eee (1) 
Ido ee } 
J=J+1 
GOle 310 
END IF 
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H=0. 00186 
TF (ETACT)=HeGis or (i) ree 
COTO] 20 
END IF 
80 DO. 90 -bSd245 ae 
WRITE (20437 %)) Yi (bh) ETA) 
20 END DO 


100 END DO 

SOE 

END 
COR KK KKK KK KK KK KR KKK KK KEK KKK KKK KKK KKK KKK KEKE KKK KKK KK KEK KKK KKK 
e SUBROUTINE DERIVative * 


CREEK KEKKEKKKEKKEKEKKEKKEKEKEKEKKEKKEKEKKKEKKEKKEKEKEKEKEKEKEKKEKEKKEKKEKEKEKKEKKKKEKKKESK 


SUBROUTINE DERIVative (ETA,N,I,Y1,Y2,YDOT, MACH) 
IMPLICIT REAL*8 (A-H,0O-Z) 
IMPLICIT REAL*8 (M) 
DIMENSION YDOT(N) , Y1(N) ,Y2(N) , ETA(N) , MACH (N) 
GOO 20210720 SOP nO a sO)) geek 
10 Y1(1I)=ETA(I) **3+(1.0D0-MACH (IT) **2) 
Ye (=a) Oly Or 2 eine) 
YDOT( 1) =iebO/¥ 251s) 
GOTO 30 
20, Yair) -ETAGt)**3 2 Ot vine Ge) 
BAO SE) AIR) ee ODO A ONIO) 
Ym@n (eh) = 0072 (an) 


30 RETURN 

END 
CO KK KK RK KK KK KK KKK KKK KK KKK KK KKK KKK KEK KK KKK KKK KEK KKEKKEKKKEKKKE 
© SUBROUTINE NEWTONCOTES * 


CREEK KEKKKKEKKE KKK KKK KEKE KKK KEKE KEKE KEKE KEKEKEKKEKEKEKEKEKKEKKKKKKKKKKKKKKKKE 


SUBROUTINE NEWTONCOTES (ETA, YDOT, Y, SAVEY, SAVED,H,N, INDEX, I) 
IMPLICIT REAL*8 (A-H,0O-Z) 
IMPLICIT REAL*8 (M) 
DIMENSION Y(N) , YDOT(N) , SAVEY(N) , SAVED (N) , ETA(N) 
GO TO (1,2,3,4,5) ,index 
1 SAVE (ly) eva) 
SAVED (I) =YDOT (I) 
ETA (I) =ETA(I)-0.25D0*H 
RETURN 
2 SAVED (I) =7.0DO*SAVED (I) +32.0D0*YDOT (I) 
ETA (I) =ETA(I)-0.25D0*H 
RETURN 
3 SAVED (I) =SAVED (I) +12.0D0*YDOT (TI) 
ETA(I) =ETA(I)-0.25D0*H 
RETURN 
4 SAVED (I) =SAVED (I) +32.0DO*YDOT (I) 
ETA (I) =ETA(I)-0.25D0*H 
RETURN 
5 Y (I) =SAVEY (I) + (H/90.0DO0) * (SAVED (I) +7.0D0*YDOT (T) ) 
RETURN 
END 
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APPENDIX B 
FORM OF a AND £6 


The expressions for a and # as introduced in [REF. 2] are 


® = +C, (1-M2)’ (iy) 


B = +C, (1-Mé2) (118) 
fem O.8=sMi<1.2 


Upon inserting egns. 117 and 118 in egns. 16 and 21 yields 


a 
a7, 
J 


Oe ae : 
yZ?% + (1-M2)* 


aria | a (120) 


a JW? + (1-MZ) 
-(1-M2) 3 


These relations are numerically integrated and plotted in 
Figs. 16 and 17. We observe the same solutions for M,21.0 as 
expected, but distinguish a different set of solutions for the 
subsonic case M,<1.0 

Applying the same procedure as in Chap. IV, we easily 
compute the expression for the boundary surfaces 


Pa 
3] MZ (y+1) 


= 4s g dx (om) 


q2o+ (1-Mé) 


oy 


The numerical integration of eqn. 121 is shown in Fig. i8@@en 
various M,, (e.g., 0.8,0.9,0.97,190,1.1,1.2). The divergmaae 
contours of the subsonic solutions appear to give 
contradictory results once investigating the pressure 


Coetiucient behavior 


= 2 

Cp = — ie go ee (SE (1-2) (122) 
Me Gy at le) 

along the boundary surface. This implies as shown in Fig. 19 


that the boundary surface within the portion of C.>1.0 feeae- 


a vealisct te eso newton 
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APPENDIX C 
FIGURES 


KSI(X) for M,=0.8-1.2 


KSI(X) 


MACH 0.9 


MACH 0.8 





( 2 3 4 5 
ee AIS 


| 


Figure 1. Numerical integration of eqns. 40 and 42. 


ETTA(Y) 





ETTA(Y) for M.=0.8-1.2 


Y - Axis 


Figure 2. Numerical integration of eqns. 41 and 43. 


E(Y) 


ECY) vs. Y for M,=0.8-1.2 
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Figure 3. Numerical solution of E(Y). 


YO vs. M, 





Figure 4. Equation 71. 


K(X) 


K(X) vs. X for M,=0.8-1.2 


: MACH 0.9 





X - AXIS 


Figu 
gure 5. Numerical integration of eqn. 72 
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SUBSONIC BOUNDARY SURFACES 
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Figure 14. Mach lines for 
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Figure 17. Numerical integration of eqn. 120. 
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Figure 19. C, for M,=0.8 . 
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